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1.0  INTRODUCTION 

This  report  documents  the  process  to  convert  right  and  left  elevation  versus  distance  course  profiles  into  triangular 
meshes  for  use  in  Road  Property  Files  (.rdf)  used  in  RecurDyn  multibody  dynamics  software.  Multiple  scripts  were  created 
in  and  processed  with  the  MathWorks  MATLAB  technical  computing  program.  Automation  was  used  to  help  decrease 
conversion  time  and  increase  productivity. 

The  road  or  terrain  on  which  a  vehicle  travels  over  must  be  defined  to  run  a  vehicle  simulation.  In  this  case,  the  simulation 
program  (RecurDyn)  required  a  specific  file  type  and  format  be  used  when  defining  the  road  or  terrain.  To  achieve  this 
format,  course  profile  data  was  required  to  be  converted  into  the  new  file  format  (.rdf). 

A  brief  overview  of  the  conversion  process  is  as  follows: 

•  Left  and  right  elevation  vs.  distance  profile  curves  were  combined  into  one  file  using  MathWorks’  MATLAB. 

•  Road  geometry  was  created  in  PTC’s  Pro/Engineer  Wildfire  5.0. 

•  Road  geometry  was  meshed  in  Altair’s  HyperMesh. 

•  The  mesh  file  was  exported  as  an  ASCII  database  file  and  imported  into  MATLAB  for  further  processing. 

•  Scripts  were  created  to  scan  the  database  file  and  extract,  format,  and  output  data  in  the  Road  Property  File 
format. 


2.0  ROAD  GEOMETRY  CREATION 

Road  Property  Files  require  node  and  element  data  in  order  to  correctly  define  the  terrain  surface  (see  Ref.  1).  The  node 
and  element  data  comes  from  data  originating  with  left  and  right  elevation  vs.  distance  profile  curves.  Combining  both  left 
and  right  curves  was  the  first  step  in  geometry  creation.  The  geometry  was  generated  in  Pro/E  to  be  used  in  Hypermesh 
to  create  nodes  and  elements  which  were  necessary  in  the  final  Road  Property  file.  The  geometry  was  created  by 
completing  the  following  steps: 


1.  A  script  file  was  created  in  MATLAB  to  combine  the  left  and  right  terrain  course  profile  files  into  one  Importing 
Blend  File  (.ibl)  The  MATLAB  script  is  included  in  Appendix  A.1.  The  elevation  vs.  distance  course  profile  files 
simply  contain  two  columns  of  data  representing  distance  and  elevation  at  that  distance,  respectively  (see  Figure 
1).  The  following  example  displays  the  left  column  as  distance  in  inches  and  the  right  column  as  elevation  in 
inches. 
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0 

-0.421332 

3 

-0.428685 

6 

-0.659071 

9 

-1.010741 

12 

-1.015889 

15 

-1.141309 

18 

-1.055657 

21 

-1.005185 

24 

-0.789697 

27 

-0.583433 

30 

-0.478471 

33 

-0.452257 

36 

-0.15139 

39 

-0.465199 

42 

-0.234502 

45 

-0.061961 

48 

-0.094372 

Figure  1  -  Sample  elevation  vs.  distance  course  profile  curve  file  containing 
points  representing  distance  in  the  left  column  and  elevation  at  that  distance 
in  the  right  column  (left)  and  a  plotted  elevation  vs.  distance  curve  (right). 


2.  The  Importing  Blend  File  contains  points  used  by  Pro/E  to  create  geometry.  The  file  defines  one  or  more  curves 
that  are  blended  together  in  Pro/E  to  generate  the  surface  geometry.  The  file  must  define  the  section  and  blend 
type  as  well  as  the  sections  and  curves  (see  Figure  2  and  Ref.  2).  In  MATLAB,  a  script  loaded  and  stored  the  left 
and  right  elevation  vs.  distance  profile  values  as  matrices.  Each  matrix  was  printed  in  its  own  section  in  the 
Importing  Blend  File  format  to  be  used  by  Pro/E  (see  Ref.  2,  Appendix  A.2,  and  Figure  3). 


Blend  File  Format 

The  imported  blend  data  file,  with  the  file  extension  .ibl,  is  in  the  following  format. 
Comments  are  contained  between  /*  and  */ 

/*  beginning  of  file  */ 

section_type  /*  The  section  type  (open  or  closed) .  */ 

blend_type  /*  The  blend  method  (arclength  or  pointwise) . 

Both  blend  types  require  the  same  number  of  curve  segments 
for  each  section.  An  arclength  blend  uses  a  general 
blending  routine  to  connect  the  sections.  The  number  of 
points  in  corresponding  curves  can  be  different  far  each 
section.  A  pointwise  blend  connects  from  point  to  point 
(point  1  in  one  curve  connected  to  point  1  in  the  other 
curve)  .  The  corresponding  curves  in  each  section  must  have 
the  same  number  of  points . */ 


begin  section 

/*  Begin  a  new  blend  section. 

This  appears  at  the  beginning  of  each  section. 

*/ 

begin  curve 

/*  Begin  a  new  curve  for  the  section. 

This  appears  at  the  beginning  of  each  curve  segment .  */ 

1  x  y  z 

2  x  y  z 

/*  The  number  is  the  point  number  (optional); 
Z  are  the  coordinate  values .  */ 

X,  Y, 

#  x  y  z 
begin  curve 

1  x  y  z  /*  The  f  irst  point  in  this  curve  equals  the  last 

2  x  y  z  point  of  the  preceding  curve.  */ 


$  x  y  z 

/»  end  of  file  */ 


Figure  2  -  Pro/ENGINEER  Blend  File  Format  help  file  describing  correct  Importing  Blend  File  format. 
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open 

pointwise 

Begin  section  !  1 

Begin  curve 

!  1 

0.000000 

0.000000 

-0.421332 

0.000000 

3.000000 

-0.428685 

0.000000 

6.000000 

-0 . 659071 

0.000000 

9.000000 

-1.010741 

0.000000 

48123.000000 

1.359609 

Figure  3  -  Sample  Importing  Blend  File  (.ibl)  defining  the 
section  type,  blend  type,  sections,  and  curves. 


3.  In  Pro/Engineer  (Pro/E),  the  Blend  from  File  tool  was  used  to  import  the  Importing  Blend  File  which  created  the 
geometry  for  the  road’s  surface.  This  tool  was  opened  by  selecting:  Insert  >  Advanced  >  Blend  from  File  > 
Surface  (see  Figure  3).  After  starting  the  tool  through  the  menu,  selecting  a  coordinate  system,  selecting  the 
Importing  Blend  File,  and  selecting  the  Material  Side,  Pro/E  will  automatically  generate  geometry  by  connecting 
the  points  in  the  Importing  Blend  File  with  a  surface  (see  Figure  4).  The  geometry  was  then  saved  in  Pro/Engineer 
as  a  Pro/E  part  file  (.prt). 


Figure  4  -  Pro/E  menu  navigation  for  the  Blend  from  File  feature  (left) 
and  resulting  road  geometry  created  by  Blend  from  File  feature  (right. 


3.0  NODE  AND  ELEMENT  CREATION 

The  Pro/E  geometry  was  created  to  automate  the  process  of  node  and  element  generation  in  HyperMesh.  The  Road 
Property  File  format  requires  that  this  node  and  element  data  be  included  in  the  file  to  define  the  road.  The  meshed 
geometry  was  exported  to  an  ASCII  mesh  database  file  (see  Ref.  3).  The  mesh  database  file  containing  node  and 
element  data  was  created  by  completing  the  following  steps: 
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1.  The  road  geometry  created  in  Pro/Engineer  was  imported  into  HyperMesh.  In  HyperMesh,  the  geometry  is 
imported  by  selecting:  Import  (green  arrow)  >  Import  Geometry  (icon)  >  File  Type:  Auto  Detect  >  Browse  for  file 
and  selecting  the  road  geometry  Pro/E  part  file  >  Import  (see  Figure  5). 


Figure  5  -  HyperMesh  menu  navigation 
for  importing  geometry. 

2.  The  automesh  tool  with  triangular  elements  was  used  to  mesh  the  surface  and  generate  the  nodes  and  elements. 
After  importing  the  road  surface  geometry,  meshing  of  the  surface  can  begin.  In  the  lower  control  panel,  select: 
2D  >  automesh  >  surfs  (select  road  geometry)  >  mesh  type:  triads  >  mesh  (see  Figure  6).  With  these  parameters 
set,  HyperMesh  will  create  nodes  and  triangular  elements  on  the  road  geometry  which  can  be  exported  to  be 
used  in  MATLAB. 
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Figure  6  -  HyperMesh  menu  navigation  for  meshing  road  surface  geometry. 


3.  The  resulting  meshed  road  was  exported  by  selecting:  Export  (red  arrow)  >  Export  FE  Model  (icon)  >  File  Type: 
HMASCII  >  Template:  HMASCII  >  File:  name  >  Export  (see  Figure  7). 
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Figure  7  -  HyperMesh  menu  navigation  for  exporting 
meshed  surface  as  HMASCII  text  file. 


4. 


The  initial  road  geometry  was  exported  as  a  HyperMesh  ASCII  mesh  database  file  in  text  format  (see  Appendix 
A. 3  and  Figure  8).  The  text  file  allowed  easy  extraction  of  the  node  and  element  data  to  be  used  in  the  Road 
Property  File. 


8a0113b55e804ca78d8822d413ddl81812ea317a790939a317e98922 

e44f 9f 822 

04d45adb23f749a2f 4e8722970f 6ab79a41badedll7355cd51fcca73 

0ab6698e0 

579ce3bf07fb3f731b821900 

END  OF  DATA 

*endgeometry ( ) 

END  GEOMETRY  DATA 

BEGIN  CONTROLCARDS 

END  CONTROLCARDS 

BEGIN  SYSTEMS 

END  SYSTEMS 

BEGIN  NODES 

*node (1,0.0 

48123.0  ,1.359609,0,0,0,1,15) 

*node (2,15.0 

48123.0  ,1.385798,0,0,0,1,7) 

*node (3,30.0 

48123.0  ,1.411978,0,0,0,1,7) 

*node(4,45.0 

48123.0  ,1.438149,0,0,0,1,7) 

*node (5,60.0 

48123.0  ,1.46431  ,0,0, 0,1, 7) 

*node (6,  75 . 0 

48123.0  ,1.490462,0,0,0,1,16) 

*node (7,75.0 

48108.92,1.496586,0,0,0,1,9) 

*node (8,75.0 

48094.83,1.479693,0,0,0,1,  9) 

*node (9,75.0 

48080.75,1.253126,0,0,0,1,  9) 

*node (10,75.0 

, 48066.67,1.100342,0,0,0,1,9) 

*node (11,75.0 

, 48052.59,0.897367,0,0,0,1,9) 

*node (12,75.0 

, 48038.51,0.612016,0,0,0,1.9) 

Figure  8  -  Meshed  road  geometry  in  HyperMesh  displaying  triangular  elements  (left)  and  a 
sample  of  the  ASCII  database  file  containing  node  and  element  data  exported  from  HyperMesh  (right). 


4.0  MATLAB  SCRIPTING 

The  technical  computing  program  MATLAB  was  used  to  create  the  final  Road  Property  file.  Several  script  files  were 
constructed  to  extract,  organize,  and  output  the  Road  Property  file  data  from  the  HyperMesh  ASCII  database  file.  The 
following  scripts  were  used: 

1 .  A  script  file  that  parsed  the  HyperMesh  ASCII  mesh  database  file,  extracting  and  storing  the  necessary  node 
information  including  node  numbers  and  the  respective  x,  y,  and  z  coordinates  (see  Appendix  A.4). 

2.  A  script  file  that  parsed  the  HyperMesh  ASCII  mesh  database  file,  extracting  and  storing  the  necessary  element 
information  including  element  numbers  and  the  first,  second,  and  third  nodes  which  comprised  each  element  (see 
Appendix  A.5). 

3.  A  script  file  that  combined  and  printed  all  necessary  information  for  the  Road  Property  file  in  the  Road  Property 
file  format  (see  Appendix  A. 6). 
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5.0  RESULTS 

Figure  9  shows  the  Road  Property  File  MATLAB  output.  Sections  of  node  and  element  data  were  removed  to  shorten  the 
text  file  for  easier  viewing. 


!  |  ROAD  DATA  FILE  | 

i 

!  Default  road  file  -  BB.rdf 

i 

METHOD 
GEN  DATA 

i 

!  Conversion  factors 

i 

X  SCALE 

i :  o 

Y  SCALE 
170 

Z  SCALE 
170 

i 

!  R  o  a  d  origin  is  located  at  the  following  global  coordinates  in  the 
!  dat  a  s  et  . 

i 

ORI  Gl  N 

0.0  0.0  0.0 

i 

UP 

0.0  1.0  0.0 

!  R  o  a  d  coordinate  system  is  oriented  with  respect  to  the  global  origin 
!by  the  following  transformation  matrix. 

i 

ORI  E  NT AT  I  ON 

1  0  0 

0  1  0 

0  0  1 

i 

! Coordi  nates  for  the  node  points  on  road 

1 

NODES 

2  0  6  8  7 

1  0.000000  48123.000000  1.359609 

2  15.000000  48123.000000  1.385798 

3  30.000000  48123.000000  1.411978 

4  45.000000  48123.000000  1.438149 

5  60.000000  48123.000000  1.464310 

6  75.000000  48123.000000  1.490462 

7  75.000000  48108.920000  1.496586 

8  75.000000  48094.830000  1.479693 

9  75.000000  48080.750000  1.253126 

10  75.000000  48066.670000  1.100342 


20687  30.022680  3037.347000  -0.848530 

ELEMENTS 
3  4  4  6  9 

1  1.000000  7.000000  6948.000000 

2  1.000000  6947.000000  6948.000000 

3  1.000000  6947.000000  6945.000000 

4  1.000000  6948.000000  6946.000000 

5  1.000000  9.000000  6940.000000 

6  1.000000  8.000000  6949.000000 

7  1.000000  6940.000000  6941.000000 

8  1.000000  6949.000000  6944.000000 


9  1.000000  6945.000000  6950.000000 

10  1.000000  6946.000000  6953.000000 

34469  1.000000  6948.000000  5.000000 


Figure  9  -  Road  Property  File  used  to  define  the  road  profile 
in  the  Redurdyn  multibody  dynamics  program 


6.0  SUMMARY/CONCLUSION 

In  order  to  simulate  a  vehicle’s  movement  over  certain  terrain  in  RecurDyn,  a  Road  Property  file  (RDF)  that  defines  the 
terrain  or  road  is  required.  The  creation  process  of  this  file  was  partially  automated  to  reduce  conversion  time  and 
increase  productivity. 

The  final  RDF  originated  from  elevation  vs.  distance  course  profile  curves.  These  curves  were  used  to  create  geometry  in 
Pro/E  that  was  meshed  in  HyperMesh  (HM)  to  create  nodes  and  elements.  Using  MATLAB,  the  node  and  element  data 
was  extracted  from  the  HM  exported  ASCII  mesh  database  file  and  printed  in  the  resulting  Road  Property  file  format.  Now, 
this  file  is  ready  to  be  used  by  the  RecurDyn  multibody  dynamics  program  in  conjunction  with  other  data  to  simulate  a 
vehicle  travelling  over  a  road  surface  or  terrain. 


7.0  APPENDIX 


A.1  -  MATLAB  script  to  create  Importing  Blend  File 


%creates  ibl  file  from  two  dist.  vs.  elevation 

curves 

%load  left  curve  and  create  x,y,z  matrix 
load  bbL.crvdet 
x  =  bbL; 

xmat  =  [ zeros ([ length (x) , 1 ]) ,  x] ; 

%load  right  curve  and  create  x,y,z  matrix,  offset  by  75  inches 
load  bbR.crvdet 
y  =  bbR; 

ymat  =  [75*ones ( [length (y) , 1] ) ,  y] ; 

%create  blendl . ibl  file 

fid  =  f open (' terrain . ibl ' ,  'w'); 

%writes  header  info  and  x  matrix 
fprintf(fid,  ' open\npointwise\nBegin  section  ! 
fprintf(fid,  ' %7 . 6f \t%7 . 6f \t%7 . 6f \n '  ,  xmat'); 

l\nBegin  curve  !  l\n'); 

%writes  header  info  and  y  matrix 
fprintf(fid,  'Begin  section  !  2\nBegin  curve  ! 
fprintf(fid,  ' %7 . 6f \t%7 . 6f \t%7 . 6f \n  ,  ymat'); 

1 \n ' )  ; 

fclose (fid) ; 

A.2  -  Resulting  Importing  Blend  File  sample 

open 

pointwise 

Begin  section  !  1 

Begin  curve  !  1 
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0.000000 

0.000000 

-0.421332 

0.000000 

3.000000 

-0.428685 

0.000000 

6.000000 

-0 . 659071 

0.000000 

9.000000 

-1.010741 

0.000000 

12.000000 

-1.015889 

0.000000 

15.000000 

-1 . 141309 

0.000000 

48123.000000 

1.359609 

Begin  section  !  2 

Begin  curve 

!  1 

75.000000 

0.000000 

-0.513719 

75.000000 

3.000000 

-0.535514 

75.000000 

6.000000 

-1.005605 

75.000000 

9.000000 

-1.388840 

75.000000 

12.000000 

-1.215158 

75.000000 

15.000000 

-1.210530 

75.000000 

18.000000 

-0 . 946456 

75.000000 

21.000000 

-0 . 859729 

75.000000 

24.000000 

-0 . 903716 

75.000000 

27.000000 

-0 . 612990 

75.000000 

48123.000000 

1.490537 

A.3  -  Sample  ASCII  mesh  database  file  (exported  from  HyperMesh) 


HYPERMESH  Input  Deck  Generated  by  HyperMesh  Version  :  10.0build60 
Generated  using  HyperMesh-HmAscii  Template  Version  :  10.0-SA1-140 
*filetype (ASCII) 

* vers ion (10 . 0build60) 


BEGIN  DATA 

BEGIN  GEOMETRY  DATA 

*begingeometry ( ) 

0000000000000000000014401f 8b080000000000000b74bc0938d45ff 43f3e635fb386b2562ada48 
2422271449911449a90845210945f 67ddf 97ac63adl445916429a2148a90648b2c65df 87313f7d7f 
5ff77efebfe7f 9bf 9f 6732f3eacceb9e73eebdaf 73ee7b065ac2ffbd58da326c9559c4818d404ba0 
5de0e4f50d3c7b7daf360facfd4c4cf877bd5209211208bbe30862ffdeb571f5blfa92b02be07f9f 
acfe201068fef7blf31f 4abbfaf 86742bc214420c8ad3ef 977fd03fe7bbladbe505535dld339acbe 
47574flfff 9722e3ff 7dfeef Odff 98b67bacfec3bbfa409474ab2ffe5df 4ab8f7fcfc53d0844ce7f 
2f 9005d3f fbe58fd4160597dfcfbcf 6dle04e0fc5fee7f f 8aa39e3fble07a512b3ca70b9ald29ceb 
8330421 60f2a0d36aba4ff870bdc84ea7f3fa207el4b8ca0fed8a01930acbe69ab07137072fc7bf 7 
bfc7bf 08d7dc597dfdf f 5cd4d5cbc46dl0cc42072ecc56ceaafclbfddf fb093cf f 8fel7f 5e3a597c 
38flef3dfea65207e2ab6755fefdd7bf71ff5dffc2fa9ffl39579ffc4bf9ffcff8947f0357e8cll9 

BEGIN  NODES 
*node (1,0.0 
*node (2,15.0 
*node (3,30.0 
*node (4,45.0 
*node (5,60.0 
*node (6,75.0 
*node (7,75.0 
*node (8,75.0 
*node (9,75.0 
*node (10,75.0 
*node (11,75.0 
*node ( 12 , 75.0 
*node (13,75.0 


48123. 
48123. 
48123. 
48123. 
48123. 
48123. 
48108  . 
48094  . 
48080  . 
,48066 
, 48052 
,48038 
, 48024 


,1. 
,1. 
,1. 
,1. 
,1. 
,1. 
92,1. 
83,  1 . 
75,  1 . 
.67,1 
.59,0 
.51,0 
.43,0 


359609,0 
385798,0 
411978,0 
438149,0 
46431  ,0 
490462,0 
496586,0 
479693,0 
253126,0 
. 100342, 
.897367, 
. 612016, 
.526117, 


,0,0 

,0,0 

,0,0 

,0,0 

,0,0 

,0,0 

,0,0 

,0,0 

,0,0 

0,0, 

0,0, 

0,0, 

0,0, 


,1,15) 

,1,7) 

,1,7) 

,1,7) 

,1,7) 

,1,16) 

,1,9) 

,1,9) 

,1,9) 

0,1,9) 

0,1,9) 

0,1,9) 

0,1,9) 


BEGIN  COMPONENTS 
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Component ( 1 , "BB_ROAD . PRT"  ,0,58,0) 
* alt surf ace (1,1,20) 

* component (2 , "autol " ,0,11,0) 

*tria3  (1, 1, 7, 6948, 6,  0) 

*tria3 (2,1,6947,6948,7,0) 

*tria3 (3, 1,6947, 6945, 694 6,0) 

*tria3 (4, 1,6948, 6946, 4,0) 

*tria3 (5, 1,9, 6940, 6949,0) 

*tria3 (6, 1, 8, 6949, 6947,0) 

*tria3 (7,1,6940,6941,6944,0) 

*tria3 (8, 1, 6949, 6944,6945,0) 

*tria3 (9, 1, 6945, 6950,6952,0) 

*tria3 (10, 1, 6946, 6953, 4, 0) 


A.4  -  MATLAB  script  to  extract  node  data  from  ASCII  file _ 

%this  script  will  parse  the  .hmascii  file  and  extract  each  node 
%number,  x,  y,  and  z  coordinates 

fid  =  fopen('BB  road . hmascii ' ,  ' rt ' ) ; 

fido  =  f open ( ' nodesall . txt ' ,  'wt'); 

%scans  each  line  of  BB_road . hmascii  and  checks  against  string:  *node ( 

%prints  any  lines  that  match  above  case  to  nodesall.txt 
while  1; 

tline  =  fgetl(fid); 
if  tline— =-l,  break,  end; 
if  strncmp (tline,  ' *node ( ' ,  6); 
fprintf (fido, ' %s\n ' ,  tline); 

end; 

disp (tline) ; 

end 

fclose (fid) ; 
fclose (fido)  ; 

%reads  nodesall.txt  and  extracts  the  node  number,  x,  y,  z  coordinates 

[node,  x,  y,  z]  =  textread (' nodesall . txt ' ,  '*node(%d  %f  %f  %f  %*[A\n]',  'delimiter'. 


^creates  nodeesformat.txt  file  to  store  formatted  nodal  information 
%prints:  total  number  of  nodes  =  length (node) 

%prints  variables:  node,  x,  y,  z 

nodesmat  =  [node,  x,  y,  z] ; 

fid2  =  f open (' nodesformat . txt ' ,  'wt'); 

fprintf ( fid2  ,  '%f\n',  length (node) ) ; 

fprintf ( fid2 ,  '%f\t  %f\t  %f\t  %f\n',  nodesmat'); 

fclose ( f id2 ) ; 


A.5  -  MATLAB  script  to  extract  element  data  from  ASCII  file 


%this  script  will  parse  the  .hmascii  file  and  extract  each  element 
%number,  node  1,  node  2,  node  3 

fid  =  fopen('BB  road . hmascii ' ,  '  rt ' ) ; 

fido  =  f open (' elemall . txt ' ,  'wt'); 

%scans  each  line  of  BB  road. hmascii  and 
%prints  any  lines  that  match  above  case 
while  1; 

tline  =  fgetl(fid); 

checks  against  string:  *tria3 ( 
to  elemall.txt 

9 


if  tline==-l,  break,  end; 
if  strncmp (tline,  ' *tria3(',  7); 
fprintf (fido, ' %s\n ' ,  tline); 

end; 

disp (tline) ; 

end 

fclose (fid) ; 
fclose (fido) ; 

%reads  elemall.txt  and  extracts  the  element  number  &  nodes  1,  2,  3 

[elem,  nl,  n2,  n3]  =  textread ( ' elemall . txt ' ,  '*tria3(%d  %f  %f  %f  %*[A\n] ',  'delimiter'. 


%creates  elemformat.txt  file  to  store  formatted  nodal  information 
%prints:  total  number  of  elements  =  length (elem) 

%prints  variables:  node,  x,  y,  z 
elemmat  =  [elem,  nl,  n2,  n3,]; 
fid2  =  f open (' elemformat . txt ' ,  'wt'); 

fprintf ( fid2 ,  '%f\n',  length (elem) ) ; 

fprintf ( fid2 ,  '%f\t  %f\t  %f\t  %f\t  %f\t  %f\n',  elemmat'); 

fclose ( f id2 ) ; 


A.6  -  MATLAB  script  that  prints  extracted  data  in  Road  Property  File  Format 


%this  script  outputs  road 

%define  conversion  factor 
xscale  =  1.0; 
y__scale  =  1.0; 
z_scale  =  1.0; 

%define  origin 
origin  x  =  0; 
origin_y  =  0; 
origin  z  =  0; 

%define  direction  for  up 
up_x  =  0 ; 
up_y  =  1; 
up_z  =  0 ; 

%define  orientation 

0  =  [1  0  0 ;  0  1  0 ;  0  0  1  ]  ; 


profile  data  in  .rdf  format 


%begin  printing  text 

fidl  =  f open (' crvdet2rdf . txt ' ,  'wt'); 


fprintf  ( fidl ,  '  \  r  !  - \n!  |  ROAD  DATA  FILE  |\n!\n!  Default  road  file 

-  BB . rdf \n ! \n ' ) ; 


fprintf (fidl,  ' METHOD\nGENDATA\n ' ) ; 

fprintf ( fidl ,  '  ! \ n !  Conversion  f actors\n ! \n ' ) ; 


fprintf (fidl , 
fprintf (fidl , 
fprintf (fidl , 


'X  SCALE\n%2 . If \n ' ,  x  scale); 
'Y  SCALE\n%2 . If \n 1  ,  y  scale); 
'Z  SCALE\n%2 . If \n ' ,  z  scale); 
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fprintf ( f idl , 

'!\n!Road  origin  is  located  at  the  following  global  coordinates 

in 

the\n!data  set 

.  \  n  !  \  n  '  )  ; 

fprintf ( f idl , 

' 0RIGIN\n%2 . If  %2.1f  %2.1f\n!\n',  origin  x,  origin  y,  origin  z); 

fprintf ( f idl , 

'UP\n%2.1f  %2.1f  %2.1f\n',  up  x,  up  y,  up  z); 

fprintf ( f idl , 

'!Road  coordinate  system  is  oriented  with  respect  to  the  global 

origin\n ! by 

the  following 

transformation  matrix . \n ! \n ') ; 

fprintf ( f idl , 

' 0RIENTATI0N\n ' )  ; 

fprintf ( f idl , 

' %d\t  %d\t  %d\n 1 ,  0)  ; 

%print  list  of 

nodes  generated  from  nod.m 

fprintf ( f idl , 

'! \n ! Coordinates  for  the  node  points  on  road\n!\n'); 

fprintf ( f idl , 

'N0DES\n%d\n ' ,  length (node) ) ; 

fprintf ( f idl , 

'%d\t  %f\t  %f\t  %f\n',  nodesmat' ) ; 

%print  list  of 

elements  generated  from  ele.m 

fprintf ( f idl , 

' \rELEMENTS\n%d\n ' ,  length (elem) ) ; 

fprintf ( f idl , 

'%d\t  %f\t  %f\t  %f\n',  elemmat'); 
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